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A Gauss-Bonnet term naturally appears in the action for gravity when one considers the existence of space 
time with dimensions more than 1+3. A variety of inflationary models can be obtained within such a framework, 
once the scale factor for the hidden dimension(s) is not constrained to be the same as that of the visible ones. In 
particular, the need for an adhoc inflaton field is eliminated. The phase space has a rich structure with different 
types of solutions, both stable and unstable. For a large class of solutions, the scale factors rapidly approach an 
asymptotic exponential form. Furthermore, sufficient inflation can be obtained for only a modest compression 
of the hidden world, if the latter is of a sufficiently large dimension. 
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I. INTRODUCTION 



The prospect of gravity propagating in a space of dimensions larger than the canonical (1 + 3) has long been an 
intriguing one. Apart from being a natural generalization of General Relativity, the idea has proffered tantalizing 
visions of both a way out of many vexing issues in high energy physics as well as that of unification of gravity with 
other forces. Starting with the efforts of Kaluza and Klein|[lll2l, passing through the formulation of theories in both 
large extra dimensions as well as small but warped extra dimensions, to gravity theories being duals to strongly coupled 
quantum field theories, it has been a riveting saga. Added to this has been the promise of String Theory as a path to 
both a consistent quantum theory of gravity as well as unification with the strong and electroweak forces. 

Naturally, this has also lead to the formulation of cosmology in space-times with a larger number of dimensions, 
and the literature abounds with a number of such scenarios. These encornpass both models wherein the visible world 
is confined to the usual (1 + 3)-dimensions, viz. the braneworld models Jj-lg], as well as those where we traverse the 
entire gamut of extra dimensions (with the latter, perforce, now being compactified with a sub-attometer radius). 

Whatever be the details of such models, it is quite apparent that the presence of such extra dimensions would make 
themselves manifest at suitably high energies, for the extra degrees of freedom associated could then be excited. Verily, 
this is the guiding principle for both collider searches as well as the study of their role in quantum corrections to low- 
energy processes. In a similar vein, their role in the dynamics of the universe as a whole would be more pronounced 
when the available energy was large, or equivalently when the distance scales were small, in other words, the era of 
the early universe. 

If we contemplate the quantum nature of gravity, it is almost self-evident that the famiUar Einstein equations cannot 
represent the whole story. Even if one started with the Einstein-Hilbert action at the classical level, quantum corrections 
would generate higher derivative terms in the effective action. The form of the latter, at first glance, is protected only 
by the symmetries of the theory, which in the present case, is general covariance. While this requirement, by itself, 
is not very restrictive, one may treat the additional terms as a power series in curvature (or, rather, in scalars formed 
out of the curvature tensor) and, in a phenomenological treatment, confine oneselves to the lowest non-trivial terms 
in this expansion. Even then, three different terms are possible. One particular combination viz. the Gauss-Bonnet 
termj?! HI Hill, is of particular interest both on account of certain interesting properties to be discussed below and for 
being the lowest order correction to the Einstein-Hilbert action as evaluated within String Theory lfl2l . In terms of the 
Riemann tensor Rf_i^(t>C,^ Ricci tensor Ri^^^ and the Ricci scalar R, it can be expressed as 

G = R^ ~ 4R^, R^" + R"""^^ ■ (1) 

Interestingly, in (1 + 3)-dimensions, this is but a topological term and adds Uttle to geometrodynamics. In a space 
of higher dimensions, this is no longer the case and the addition of this term to the action can lead to non-trivial 
modifications to the equations of motion, thereby giving rise to interesting effects. In this paper, we will analyze some 
of the features of such a theory in context of the early universe. 

In section 2, we look at the evolution equations and the conservation equation, following by section 3 which lists 
all the cosmological solutions and discusses the behaviour of the scale factors for both the visible and the extra 
dimensions. In the end, i.e., section 4, we summarize the paper and give conclusions. 



II. EVOLUTION EQUATIONS 



The Einstein equations, as derived for the Einstein-Hilbert action are second-order in the derivatives of the dynam- 
ical variable, viz. the metric. This is true for any dimension of spacetime. Since the dynamics of the Universe today 
is well explained by this theory, it is paramount that any deviation from the Einstein-Hilbert action should essentially 
disappear in the context of the present-day universe atleast to first order In particular, this motivates the thought 
that even for higher dimensional theories, the equations of motion need to be second ordefl It restricts the form of 
the action and the simplest term that one can write beside the Einstein-Hilbert and cosmological constant term is the 
Gauss-Bonnet term. Such a deformation has the further advantage in the cosmological context that on account of it 
being a quadratic in the curvature, it could have played a significant role in the early history of the universe, while its 
effect today would be negligible. 

We consider the higher dimensional Gauss-Bonnet actionlflsi fTill as 



4+D 



X\/—g 



^rn ^(Ai? +eG) 



(2) 



where R and g correspond to the {1 + 3 + D) dimensional curvature scalar and determinant of the metric tensor 
respectively. Although we have not added a cosmological constant term explicitly, we can always if intended, do so 



While this requirement is not a strict one, it is certainly an useful one, and also serves to eliminate ghosts, which are generic in higher dimensional 
theories. 



3 



by adding an appropriate piece in the 'matter' lagrangian £. M defines the scale of quantum gravity in the entire bulk, 
and is related to A/pianck (defined in 4-dimensions) through VdM^^"^ — M^^^^^^, where Vd is volume of the extra- 
dimensional subspace. One would expect that Vd ^ M^^ and, thus AI < A/pianck- The quantity e parametrizes the 
weight of the Gauss-Bonnet term and can be reexpressed as e = ^M^"^ where 7 is a dimensionless constant. If the 
Gauss-Bonnet term is the result of quadratic corrections, then one would typically expect I7I < 1. The two choices 
A = +1(— 1) correspond to attractive(repulsive) gravity in the Newtonian Umit. While A = — 1 is not unreasonable in 
an epoch where the Gauss-Bonnet term dominates, clearly it is untenable at late times (i.e., the present epoch) where 
the curvature is small and the standard Einstein-Hilbert action needs to be recovered. We shall, then, assume A = +1. 
It should be noted that String Theory predicts that 7 < 0, a resultigl [l2l that we shall adopt. 

What we can demand purely from observations is the homogeneity and isotropy of the normal 3— dimensional space. 
If this is the only constraint to be satisfied, the line elment can be expressed as 



dr^ 



1 - fcir2 



r dil + 2goidtdx + gjxdx dx 



(3) 



We will, however, resort to a simplification in which the hidden space is also homogeneous and isotropic. Neither 
the observed homogeneity and isotropy nor any other conclusion about observable cosmology that we shall draw 
is contingent on this assumption. On the other hand, this serves to constrain the degrees of freedom in the theory, 
thereby enhancing its falsifiability. The scale factor a{t) will, then, be a function only of time. Similarly, the extra 
dimensional sector of the line element can also be expressed in a form similar to the normal dimension with its own 
scale factor, b{t). While both these sub-spaces are isotropic among themselves, the overall space-time is manifestly 
anisotropic. This feature is built into the formalism by allowing a{t) and b{t) to be different. Such a space-time which 
has one temporal dimension, 3 'visible' spatial dimensions and D extra spatial dimensions, can be described by the 
line element, 



-df + a^it) 



dr^ 



1 - kir- 



r^dn 



b\t) 



dR^ 



1 - /C2i?2 



R^dVlD^ 



(4) 



We denote the time dimension with super(sub)script '0'. Whereas we reserve lower-case Roman indices (i,j = 
1, 2, 3) for the normal spatial dimensions, upper case Roman indices denote the extra dimensions and take values 
J, J = 4, 5, . . . , D + 3. The observable universe is well described by a vanishing spatial flatness and hence ki = 0. 
For the sake of simplicity, we likewise assume the extra dimensional curvature to be vanishing i.e. ^2 = 0. In other 
words, our spacetime manifold is described by Ad^ x T^. 

Starting from equation (|2]i, the Euler-Lagrange equation for the metric would read 



1 



ii/2+D 



(5) 



where is the 1 + 3 + Z? dimensional stress-energy tensor for the matter content, and the generalised Einstein tensor 
G^^ incorporates the contribution from the Gauss-Bonnet term. Its non-zero components are given by 
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While the Gauss-Bonnet term gives rise to terms quadrilinear in the scale factors (a(t) and b{t)) or their derivatives, 
they still contain only upto the second derivatives. (This feature is not particular to the cosmological solution, but is 



- Note that k2 ^ 0, for example, A/'* X would, typlically, require non-trivial energy-momentum tensor. 
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also retained in the general case.) It is worth noting that, for D = 1, eq.® is independent of b{t) while this is not true 
for D > 1. This distinction is of profound importance and leads to some very interesting dynamical features that some 
earlier efforts ifioll . concentrated as they were on D = 1, had missed. Indeed, the addition of the sixth, seventh and 
eigth dimensions, each brings into play a new term in the dynamics. Adding further to D does not introduce any new 
term; instead, it would alter both the overall weight of the Gauss-Bonnet term relative to the Einstein-Hilbert term, and 
also subtly shift the balance between the various 0{e) terms in the equations above, thereby affecting the dynamics. 

Note that our assumption of a metric of the form of equation ©, or in other words, homogeneity and isotropy of 
the universe separately in the visible and the extra dimensional subspaces require that the stress-energy tensor be of 
the form 

Ti; = diag{-p, Pa, Pa, Pa, Pb, Pb) (9) 

Here, p is the energy density of the fluid which is present everywhere. However, the pressure it exerts on the two 
subspaces are different in general. 

Conservation of energy-momentum (T^^ = 0) implies 

^ipa'b^) + Pab^^ + Pba'^ = (10) 
at at at 

For the sake of simplicity, the equations of state are assumed to yield Pa = WaP and = w^p, where Wa and Wb are 
constants. On solving equation ( fTOl i. we, then, have 

3(l+«)„) ,D(l + u,,) 

Here, po is the energy density corresponding to the epoch when a = ag and b — bo. While we could have generalized 
equation ^ to include multiple fluids, this would have added Uttle to our analysis. Hence, we desist from such 
complication. Indeed, as we shall see later, the matter content has relatively little importance. 



III. COSMOLOGICAL SOLUTIONS 



A. General features 



The evolution of the Universe is governed by the Einstein's equations (with the conservation equations being built 
into it.) As can be seen, we have five variables and six equations, namely the three Einstein equations (161718b . the 
conservation equation (fTOl i and the two equations of states. However, once we Eq. [TO]is considered, only one out of 
equations ( I7l8l l is independent. 

We start the evolution att = M^^. While the classical equations are involved for t < M^^, the particular choice 
is just one of the convenient choices and does not have any influence on the dynamics as a result of time translation 
invariance. It will be convenient to express all the quantities in terms of dimensionless units, namely 



= - = ^ ^ 77 = 777 

a aAl b bM 



P p = Pa P = 

— ^.^j_n J^b — 



b' b 

(12) 



Here, prime represents derivative with repect to t. The initial value of t is denoted by tq. Note that for fci 2 = 0, the 
equations are independent of the absolute values of a{t) and b{t) but depend only on the values relative to those at 
present time. This is reflected by the fact that a{t) and b{t) enter into the equations(for fci.2 = 0) only as logarithmic 
derivatives. 

We can re-express the Equations (I6l7l8t in terms of dimensionless quantities as 



-p = -X[3Hl + 3D Ha Hb + °C2 H^] +12^ [D HI Ht + 3 ^Ca Hi + 3 ^Cg Ha Hi + H^] (13) 

Pa - H'^ h {Ha ,Hb)+HU2 {Ha , Hh) + h{Ha, Hb) (14) 

n = H' UHa, Hb) + H' fJHa, Hb) + hiHa, Hb) (15) 
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where 

fi{Ha,m) = 2 {-X + j[ADHaHi, + i''C2H^]) 

f2iHa,Hb) = -DX + 4-/ [DH^ + S^'C^Hi +A''C2HaHb] 
h{Ha,Ht,) = -X [3 + D+^C^Hi + 2DHa Ht] 

+ -f[8DH^Hb-2Di3- 5D) Hi Hi + i{D-l) ^Cg H^ + 8 D Ha H^] 

(16) 

h{Ha,Hb) = -{D - 1) X + 12^ [""-^Cs Hi +2''-^C2 Ha Hb + iD-1) Hi] 
h{Ha, Hb) = -3 A + 12 7 [Hi + 2 {D - 1) Ha Hi, + Hi] 
h{Ha,Hi,) = -X [6Hl + ''C2Hl + 3{D-l)HaHi,] 

+ 12 7 [H^ + H^ + {D-1) {3 Hi H^ + {2D - 3) Hi Hi + ^-^Ca Ha H^}] 

Note that eqn.(fT3]l is invariant under the joint transformation (Ha. Hh) {—Ha,—Hb). Similarly, /i(_ffa, ^6 ) ~ 
fi{—Ha, —Hb). This is reminiscent of the analogous property of the (1 + 3)-dimensional case. Thus, for an empty 
universe, both expanding (Ha > 0) and contracting {Ha < 0) solutions for the scale factor of the visible sector are 
equally possible and the growth (or otherwise) of the universe is determined by the initial conditions. We will return 
to this point later 

It is instructive to examine the structure of eqn.(fTjt alone, keeping dynamics aside for the moment. Being cubic 
in Ha, this leads to three solutions for a given (7, p, Hb) combinations. In other words, there are three branches for 
the solutiorQ. In FiglT] we depict these branches as plots of Ha{Hb) for specific values of 7 and p. While p = 
corresponds to an empty universe, p — 1 would indicate a cosmological constant comparable to the scale of quantum 
gravity, or in other words, the limit of applicability of the theory. Thus, during its evolution, the energy density of 
the universe must lie somewhere in between these limits. Depending on the values of (7, p, Hb), the solutions for Ha 
could either be all real, or one real and a complex-conjugate pair Complex solutions are, of course, not admissible on 
physical grounds, and this is indicated in FiglUby the fact of there being only one solution for certain values of Hb. 
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FIG. 1. Curves correspond to different 7 values for p — 0.0[Left Panel] and p = 1.0[Right Panel]. 

Irrespective of the nature of the dynamics, it is obvious that, for a given p{t), the pair {Hb, Ha) must lie on a curve 
that is bounded by the curves for p = and p = 1. The extent to which they can move away from a specific curve 
is governed by the evolution of p{t). As Fig[T] shows, this dependence on p can be significant, especially for small 
{Ha,Hb). However, if for some reasons, such a configuration is reached only at late times, this dependence would 
turn out to be irrelevant. Even a large primordial p would be seen to be diluted to insignificant amounts. 

Before we solve for the dynamics, it is useful to dwell upon a few more issues: 

• Note that we can never have both iJ^ 5 > for p = 0. This is a consequence of the signs of 7 and A, dicated 
as they are by other physical considerations. This, though is not true for p > 0. However, since p quickly 
decreases in magnitude to almost a vanishing value, if the visible universe is expanding at a given epoch, the 
hidden dimensions are contracting and vice versa. 



^ While it may seem that with eqn.(T3) being quartic in Ht, there should four branches instead, clearly such an argument is faulty. Note that, for 
D = 2, there are only two branches 
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• The above result has a striking consequence. Unless the universe started with a configuration opposite to what 
we have today — namely a large &(t = 1) and a small a(T — 1) — we must have had a relatively small \Hii{t)\ 
during the entire history of the universe. In particular, the inflationary epoch should have had Ha ^ \Hb\- This 
immediately tells us which of the branches in FiglUshould the initial reside on. 

• It is easy to appreciate that a large value for \Hi, \ in the present epoch would be inadmissible as this would lead 
to the effective cosmological parameters changing rapidly. Moreover, since the theory is applicable only for 
6(r) > 1, the contraction of the hidden dimensions must have essentially stopped at some distant pasQ. Since 
Hh w implies that i?^ = p/3 A (a result that one would have obtained in standard (1 + 3)-dimensions as well), 
we must have A > as we have already discussed. 

Since it is possible to choose initial conditions such that |i?6(r = 1)| <C Ha{T = 1), it is now clear (courtesy FigJT) 
that this hierarchy would be maintained. This has the immediate consequence that /5(t) would be a fast decreasing 
function of r (see eqnfTTTi as long as Wa > —1 and Wb is not very largeQ Since this inequahty is satisfied by all 
ordinary matter, it is quite apparent that the universe quickly settles down to a phase where the matter density (or 
the pressure thereof) has little bearing on the evolution of scale factors. This resullQ helps us understand better the 
assymptotic nature of the solutions to the dynamical equations that we now turn to, in particular, the lack of sensitivity 
to the matter content. 
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FIG. 2. Behaviour of the Ha and Ht with the rescaled time r for 7 = —0.1 [left panel] and 7 = —1.0 [right panel]. In either 
case, we have chosen the parameters Wa = 1.0/3.0, wt — 0.0, i.e. a pure radiation confined to (1 + 3)-dimensions, and the initial 
conditions p(r = 1) = 1, Ht{T = 1) = —1.0. 



For a Unear equation of state (as we have assumed), we can easily eliminate p and Pa by combining eqs.( fT3^T4l i. 
With Ha = Ha{Hb), what remains to be done is to solve a single first-order non-linear differential equatioj. In 
Fig[3] we illustrate the Hubble parameters for certain specific choices of the constants 7, D and the initial condition 
Hb{To) ~ —1. In each case the particular branch of the solutions Ha{Hb) has been chosen so as to lead to acceptable 
amount of inflation. Note that, in each case, Hb initially increases from Hb{To) = —1, only to reach a fixed point. 
Courtesy eqn.(fT3Tl, a similar fate befalls Ha too, with it decreasing to a positive fixed point value. Increasing D has the 
twin effect of increasing \H^^"'^p\ (decreasing \H^''^'^^^\), while delaying the epoch of settling at the fixed point. The 
dependence on I7I is analogous to that on D. This can be appreciated by realizing that the effect of the Gauss-Bonnet 
term becomes more pronounced with increase in D and decrease in |7|. 

It is useful to consider, here, the senstivity of the evolution to the parameters Wa.b as well as the initial condition 
p(to). Let us consider the latter point first. As we have argued above, a comparison of the two panels of FiglT] 
demonstrates that there is little difference between p = 1 and p = 0. On the other hand, naturalness demands that 
jO ^ 1. Thus, even though changing p(to) does have a bearing on the evolution, the effect is discernible only for a 
small time scale (the transient region) and quickly becomes insignificant. For precisely the same reason, the effect of 
Wa,b is marginal too, as long as that the energy density does not remain unchanged from a significant starting value, 
(i.e., unless u>a,& conspire so that p{t) ^ 0{l) over a significant range r ^ tq). 

A constant positive Ha (negative Hb) implies exponentially inflating visible universe (deflating hidden dimensions). 
In Figl3j we illustrate the resultant scale factors. As mentioned earlier, all evolution is independent of the initial value 



This conclusion can be evaded only if the universe had started with an immensely large 6(7"). 
' In other words, the contraction in 6(r) should not offset the expansion in a(r) in the evolution of p{t). 
^ Note that we have not assumed this simplification in our actual calculations. 

^ One could also have solved eqns. )14l & |15) for H'^ ^ and, thereby obtained not only dHa/dHi,, but also Ha{t) and Hi,{t). The two procedures 
are eniiivale.nt. 
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of the scale factors. Understandably, at late times a(r) continues to increase with the time r, while &(t) continues to 
decreas^ Note, though, that a(r) cannot be inflating forever Nor can 6(t) decrease continually, for that would either 
imply a current value of ^(t) to be smaller than the scale of quantum gravity (a regime beyond the validity of our 
treatment) or a very large initial value 6(to), which, though allowed in principle, is aesthetically unpleasing. Thus the 
goal should be a large value of the ratio i?cnd = a(''cnd)/&(7'cnd), where Tend denotes the epoch of the end of inflation. 
Without going into the mechanism of ending inflation, it is clear from Figj3]that i?cnd increases monotonicall>0 with 
D, and, thus, a larger D proffers, in some sense, a better solution. 




FIG. 3. Behaviour of the normalized scale factors a(r) and 6(r) with the reseated time r for _ff(,(1.0) = —1.0, Wa = 1.0/3.0 and 
Wb — 0.0. The left and right panels correspond to 7 = —0.1 and 7 = —1.0 respectively. 



B. Asymptotic dynamics 

The numerical solutions obtained in the preceding section suggest that, after a transient behaviour, the evolution, at 
late times asymptotes to exponential growth (decay) for the visible (hidden) sector We examine, now, whether this 
feature is a generic one or is it specific to a class of initial conditions and/or parameters. 

The existence of exactly exponential solutions for both a(T) and &(t) would, of course, be contingent on having 
H'^ = H'f^ = 0. Now, H'^ J can be algebraically solved for (in terms of Ha.b) using eqns.(fl"4MT5ll. Vanishing derivatives 
of bubble parameters, thus, translate to curves in the Ha-Hi, plane (each corresponding to specific values of Pa and 
Pb). The overlap of these curves with those corresponding to the constraint equation (fTsT l — see, for example, FiglT] — 
albeit for an allowed combination (^p{t), Pair), Pb{T)^ would, then, indicate the fixed points. While this entails 
solving the full set of the dynamical equations, it is instructive to consider a special choice, namely a constant p{t). 
For the linear equation of state that we have assumed, this essentially implies a cosmological constant- In Figlll we 
illustrate these fixed points for a particular configuration, namely (Z) = 4, 7 = —1.0). 

Not all the fixed points would be stable. Stability, in this context, implies that time evolution from the neighbourhood 
of this point should bring Ha^b to this point rather than take it away from it. In other words, for a small displacement 
from the fixed point, the derivatives H'^ ^ should point towards the point. In FigHl small triangles (circles) denote the 
stable (unstable) fixed points. Note that each stable point {Ha,Sb) is paired with an unstable one at {—Ha, ^Hb)- 
This pairing can be understood by recalling that fi{Ha, Hb) = fi{—Ha, —Hb) (see eqnfTSll. and, hence the derivatives 
must satisfy H'^ bii^a-b}) — ~-^a bii'^^a.b})- Once a stable fixed point {Ha = a,Hb = (3) has been reached, the 
system would continue in that state, leading to exponential growth (decay) thereon. 

At this stage, it is imperative to distinguish the case of the cosmological constant from the more general scenario 
of an evolving p{t). Assuming a non-trivial fixed point 7^ (0,0) exists, the former must satisfy eqn.(fTjt. 

On the other hand, the matter energy density (other than that for a cosmological constant) changes with time. This, 
immediately, demonstrates that no such fixed point can exist in the generic case. Nonetheless, as we argued earlier, for 
Wa > —1, the matter density quickly falls to insignificant levels (as long as we do not have Wb ^ —1), and thus, for 
such matter content, the asymptotic behaviour would be almost indistinguishable from that due to an empty universe. 



* We postpone the discussion of monotonicity at all times r until later. 
' This could have been gleaned from Fig|2]as well. 

It is interesting to note the possibility of having a solution with Wa,6 7^ 1. but with the Hubble parameters satisfying a'^ (i+'^n) (i+iub) ~ 

constant, at least asymptotically. 
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FIG. 4. Phase space plot of Ha and Hi, for D = 4. Curves in this figure correspond to p = O.0(solid line) and p — 1.0(dotted line) 
for7 = -1.0. 
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0.31 


-0.37 


0.83 


-0.31 


-0.37 




-0.1 


2.9 


-0.25 


-0.09 


-2.9 


0.25 


-0.09 






-1.85 


0.69 


-0.37 


1.85 


-0.69 


-0.37 



TABLE I. Different analytic solutions for p = 



Reverting to the cosmological constant, it is clear that, for D = 3, the symmetry between the visible and the hidden 
sectors guarantees that every (un)stable fixed point {a, /3) would be accompanied by another of the same ilk, namely 
a). For I? > 4, on the other hand, no such symmetry is present, and we would, in general have two rather different 
stable points (each with its unstable mirror image). For D = 2, the aforementioned three solution branches collapse 
to two, and the two stable fixed points merge. These features have been exhibited in TableU] 

For a given D, if 7 is changed, the asymptotic values (a, /3) do change, but the ratio 5 = fi/^ remains unchanged 
(see TableU). This can be understood by noting that, for an empty universe, eqns. (fT3] - [T5] l. can be trivially recast in 
terms of the scaled variables J—^ Ha and J—jHt,. The very same property dictates that the values a and B are 
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proportional to 1 / v 1 7 1 ■ If we denote the epoch of reaching the fixed point by t/ , then, at a later time r > r/ , we have 

a(r) _ a{Tf) 



b{r) b{Tf) 



exp[(«-/3)(T-T/)] (17) 



Thus, to maximize the growth of the ratio of the scale factors, we need to maximize the difference (a — (3). Clearly, 
this is aided by decreasing |7|. While this may seem paradoxical at first, this is not so given the fact that the limit 
7 ^ is not straightforward and has to be taken with car This is demonstrated by Table |I] 

The dependence of {a, (3) on the number of hidden dimensions D is more subtle. Analytically, it is best understood 
in the limit of very large D. Approximating, say, ^C„ « D"^/nl etc., the pure-gravity version of eqns. (IT3l - IT5t can 
now be recast in terms of [Ha, Hi, = D Hi,) with all dependence on D being absorbed in Hi,. This immediately 
implies that f3 should scale as D^^ while a should be free of any Z3-dependence. Of course, such a scaling is exact 
only in the infinite-D limit, and corrections to this behaviour should be apparent for moderate D values. A perusal 
of TableUshow, nonetheless, that f3 cx D^^ is a fairly good approximation, while the relative variation in a is much 
smaller Indeed, even these residual dependences can be understood analytically if a 1 /D expansion is performed. For 
our purposes, though, such an exercise is of limited use and we shall desist from it at the current juncture. Once again, 
considerations analogous to those of eqn.(fT7]i tells us that a larger D value is more suitable for effecting large inflation 
with a moderate contraction of the hidden dimensions. Once again, TableUclearly demonstrates this. 



C. End of Inflation 



One issue that we have not addressed so far is the mechanism to end inflation. In different models of extra- 
dimensions motivated inflationary scenarios iflsl - flQll . the issue of the end of inflation and reheating is addressed in 
different ways. As the scale factor of the extra-dimension decreases, it is expected that a stage is reached when quan- 
tum gravity effects become important. Our analysis, however, is of classical nature, and, hence, is not equipped to 
address this issue directly. Several speculations can be made though. For one, it is conceivable that quantum gravity 
itself would generate a pressure that would arrest such an interminable contraction, thereby forcing a{t) to come out 
of the inflationary phase to a FRW-like one. 

While this may seem to be a dissatisfactory solution, let us consider if such a situation can be reached dynamically 
with some (relatively) well-understood physics input. If b{t) is to reach a constant value, we must impose Hi, = and 
i?^ = in eqns. (113114115b . A consistent solution needs 

-_ 3 - , y/g - 967Pb + 727P,(2P,'^T3) 

p--2^a+ ^ • (18) 

Of course, such an equation of state seems extremely unnatural and contrived. Moreover, assuming this to be an 
identity would not admit the dynamics that we have seen so far, and, thus, would negate the very reason for our 
analysis. Rather, what we need is a equation of state that, at some suitably late epoch, smoothly flows to the above 
or at least to a close approximation thereoQ Even this flow seems a bit contrived if one considers only a single 
matter component. However, on the inclusion of different fields, with differing kinetic terms and potentials, such an 
eventuality in the evolution of the total stress-energy tensor can be managed with relative ease. 

Yet another possibility arises if we treat the coefficient of the Gauss-Bonnet term (e) not as a constant but a dynamical 
quantity. This is not an unreasonable assumption for the Gauss-Bonnet term and is supposed to have arisen as the result 
of quantum corrections. Once e is assumed to be a field, it will, naturally, evolve with time. If the potential for e be 
such that e = is a (local) minimum, the evolution equations would naturally force the visible universe to come out of 
the inflating phase to a decelerating era. This issue and other associated questions such as reheating and the generation 
of density perturbations will be addressed in a subsequent work. 



IV. CONCLUSIONS AND DISCUSSIONS 



In the context of a higher dimensional cosmology, we have, in this paper, given a mechanism to produce inflation 
using the Gauss-Bonnet term which provides a natural extension of the Einstein-Hilbert action. Unlike several common 
models of inflation where a field (called the inflaton field) is required for an inflationary phase ||20| - |23| in the Early 
Universe, in our case, we have eliminated the necessity of such a field. The universe goes into an inflationary phase 
just from the interplay of the extra dimensions and the Gauss-Bonnet term in the action. 

The parameter space characterized by the the expansion rate in the normal and extra dimensions. Ha and Hi, , 
exhibits a rich structure with a variety of possible solutions depending on the initial conditions. There are both stable 



Note that 7 (or, equivalently e) essentially introduces a new scale in the theory. For example, the aforementioned scaling of i, could, instead, 

be understood in ternis of rescaling time by the same factor of y'— 7- 

The approximation has to be good enough so that 6(t) does not change appreciably. 
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and unstable fixed points. If, for simplicity, linear equations of state are assumed, then the matter density dies down 
fast with time, as long as Wa > —1 and Wb < 1- Thus, asymptotically, the system reaches the limit of vacuum 
solution, which in our case turns out to be exponential behaviour of scale factors with time. More importantly, if one 
of the scale factors increase with time, the other decreases. This is a very desirable situation as we would like the 
normal dimensions to expand while the extra dimension should become smaller A large expansion is achievable for 
a significantly wide range of the dimensionless Gauss-Bonnet parameter 7. Indeed, a 7 small enough to be consistent 
with a quantum-mechanical origin is admissible on this count. And while the asymptotic values of Ha and Hi, do 
depend on the Gauss-Bonnet parameter, their ratio is essentially insensitive to it. For a given value of 7, the expansion 
rate of the visible dimensions (Ha) depends only weakly on the number of extra dimensions, D. On the other hand, 
Hb has a stronger dependence on D. In fact, for large D, Hf, is roughly inversely proportional to D^^. 

The end of inflation, the generation of density perturbations and reheating are some of the open questions in this 
context. While we have proferred some possible solutions to the first of these, a detailed mechanism fulfilling all 
constraints needs to be worked out. 
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